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I. INTRODUCTION 
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We investigate entanglement production in a class of quantum baker's maps. The dynamics of 
these maps is constructed using strings of qubits, providing a natural tensor-product structure for 
application of various entanglement measures. We find that, in general, the quantum baker's maps 
are good at generating entanglement, producing multipartite entanglement amongst the qubits close 
to that expected in random states. We investigate the evolution of several entanglement measures: 
the subsystem linear entropy, the concurrence to characterize entanglement between pairs of qubits, 
and two proposals for a measure of multipartite entanglement. Also derived are some new analytical 
formulae describing the levels of entanglement expected in random pure states. 
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\ The introduction of "toy" mappings that demonstrate essential features of nonlinear dynamics has led to many 
■ insights in the field of classical chaos. A well known example is the so-called baker's transformation pj. It maps 
the unit square, which can be thought of as a toroidal phase space, onto itself in an area-preserving way. Interest in 
the baker's map stems from its straightforward characterization in terms of a Bernoulli shift on the binary sequence 
that specifies a point in the unit square. It seems natural to consider a quantum version of the baker's map for 
, the investigation of quantum chaos. There is, however, no unique procedure for quantizing a classical map; hence, 
f^S different quantum maps can lead to the same classical baker's transformation. 

I/"") \ Balazs and Voros |2J were the first to formulate a quantum version of the baker's map. This was done with the help 
of the discrete quantum Fourier transform. Subsequently, improvements to the Balazs- Voros quantization were made 
by Saraceno Q , an optical analogy was found Q , a canonical quantization was devised 0, 0| , and quantum computing 
realizations were proposed A related quantum baker's mapping on the sphere has also been defined 0. 

More recently, an entire class of qua ntum baker's maps, based on the 2 Ar -dimensional Hilbert space of N qubits, was 
Qh, proposed by Schack and one of us [jjj] . The qubit structure provides a connection to the binary representation of the 
' classical baker's map. This connection comes through the use of partial Fourier transforms, which are used to define 
orthonormal basis states that are localized on the unit-square phase space. The n-th partial Fourier transform, which 
2 ■ acts on N — n of the qubits, defines orthonormal states that are localized at a lattice of phase-space points specified by 
qh| n position bits and N — n momentum bits. Each state is localized strictly within a position width 1/2" and roughly 
within a momentum width l/2 N ~ n . The n-th quantum baker's map in the class takes the states defined by the n-th 
partial Fourier transform to the states defined by the (n — l)-th partial Fourier transform. This action decreases the 
number of position bits by one, while increasing the number of momentum bits by one, thus mimicking the stretching 
and squeezing of the classical baker's map. By using this procedure, one obtains N different quantum baker's maps, 
one map for each number of initial position bits (or initial momentum bits). The Balazs- Voros quantization is but one 
member of this class, corresponding to a single initial position bit (n = 1). The map at the other extreme (n = V), 
which has no initial momentum bits, is easily shown to be unentangling. 

The classical limit of this class of baker's maps has been investigated by two groups 0, 0]. Tracy and one of 
us |l2fl found that if the number of initial momentum bits is allowed to approach infinity as the overall number of 
qubits goes to infinity, the classical baker's map is recovered. This result is consistent with the findings of Soklakov 
and Schack pd| . In contrast, if the number of momentum bits is held constant as the number of qubits increases to 
infinity, a stochastic variation of the classical baker's map is created The simplest such limit, which holds the 
number of initial momentum bits constant at zero (n = N), follows a sequence of completely unentangling quantum 
baker's maps. 

Our curiosity now poses the following two questions. What is the entangling power of all the quantum baker's 
maps? And what role, if any, does entanglement play in the the classical limit? This paper focuses, for the most 
part, on the first of these questions, investigating the entangling power of the Schack-Caves class of quantum baker's 
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maps. Previous investigations of entanglement in quantized chaotic systems, for the mos t part, have dealt with the 
correlations induced by coupling two or more independent systems together |liilMliiillflll3lialia|2^ Our 
approach here is quite different: each of our quantum baker's maps lives in a Hilbert space with a qubit tensor-product 
structure; strings of qubits form a natural basis, anchoring Hilbert space to the corresponding classical phase space, 
and the quantum dynamics of our baker's maps is defined explicitly in terms of this connection. We therefore expect 
an intimate relationship between the dynamics and the multipartite entanglement induced amongst the qubits. Using 
different approaches, the dynamics of entanglement in qubit bases was recently investigated in |23l l24| . 

In order to calibrate the entanglement produced by the quantum baker's map, we compare it with the entanglement 
of random pure states drawn from the appropriate Hilbert space. Thus a by-product of our investigation is to 
derive some new exact formulae describing the levels of entanglement expected in random pure states. As measures 
of entanglement, we examine in detail the subsystem linear entropy, deriving formulae for the variance and third 
cumulant, and two proposals for a multipartite entanglement measure, where formulae for the mean and variance 
are given. Pairwise (mixed-state) entanglement between two qubits drawn from N qubits is investigated numerically, 
using the concurrence as the entanglement measure. 

The paper is organized as follows. In Sec. [H] we introduce the baker's map, both in classical and quantal form. 
Section IlIII is devoted to discussing the measures of entanglement and the entanglement of typical pure states. In 
Sec. lIVI we explore the entangling power of the quantum baker's maps. Finally, in Sec.[V] we provide a brief discussion 
of our results. 



II. THE QUANTUM BAKER'S MAP 

The classical baker's map is a standard example of chaotic dynamics 1]. It is a symplectic map of the unit square 
onto itself defined by 

q n +i = 2q n - [2q n \ , (1) 

Pn+i = ( Pn + [2q n \) /2 , (2) 

where q,p € [0, 1), [^J is the integer part of x, and n denotes the n-th iteration of the map. Geometrically, the map 
stretches the unit square by a factor of two in the q direction, squeezes by a factor of a half in the p direction, and 
then stacks the right half onto the left. 

Interest in the baker's map is due mainly to the simplicity of its symbolic dynamics. If each point of the unit square 
is identified through its binary representation, q = 0-SiS2 • • • = J2T=i s k2~ k and p = 0-SoS-i . . . = J2T=o s -*2 -1 
(si G {0, 1}), with a bi-infmite symbolic string 

s = ... s_ 2 s_iSq • sis 2 s 3 . . . , (3) 

then the action of the baker's map is to shift the position of the dot by one point to the right, 

s -> s' = . . . s_ 2 s-is si • s 2 s 3 .... (4) 

For a quantum-mechanical version of the map, we work in a D-dimensional Hilbert space, TLd, spanned by either the 
position states \qj), with eigenvalues qj = (j + 0)/D, or the momentum states \pk), with eigenvalues pk = {k + ot)/D 
(j, k = 0, . . . , D— 1). The constants a,(3 S [0, 1) determine the periodicity of the space: | = e~ 2 ™ a \qj), \pr+d) = 
e 2m P\pk). Such double periodicity identifies TLd with a toroidal phase space. Periodic boundary conditions correspond 
to a = /9 = 0, and anti-periodic boundary conditions to a = (3 = 1/2; because of other symmetry considerations, 
these are the only two cases of interest. The vectors of each basis are orthonormal, (qj\qk) = (Pj\Pk) = Sjk, and the 
two bases are related via the finite Fourier transform, 

( qj \F D \q k ) = ( qj \p k ) = ^=e^ h . (5) 

For consistency of units, we must have 2irhD = 1. 

The first work on a quantum baker's map was done by Balazs and Voros 0. Assuming an even-dimensional Hilbert 
space with periodic boundary conditions, they defined a quantum baker's map in terms of the unitary operator B 
that executes a single iteration of the map. To define the Balazs- Voros unitary operator in our notation, imagine 
that the even-dimensional Hilbert space is a tensor product of a qubit space and the space of a (_D/2)-dimensional 
system. Writing j = x(D/2) + j', x S {0,1}, we can write the position eigenstates as \qj) = \x) ® \j'), where the 
states \x) make up the standard basis for the qubit, and the states \ f) are a basis for the (D/2)-dimensional system. 
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The state of the qubit thus determines whether the position eigenstate lies in the left or right half of the unit square. 
The Balazs-Voros quantum baker's map is defined by 



B = F n o 1 



(i 2 ®F D ; 2 ), (6) 



where I2 is the unit operator for the qubit, and F D / 2 is the finite Fourier transform on the (D/2)-dimensional system. 

The unitary B does separate inverse Fourier transforms on the left and right halves of the unit square, followed by 
a full Fourier transform. Later Saraceno [j| improved certain symmetry characteristics of this quantum baker's map 
by using anti-periodic boundary conditions. 

Taking again the anti-periodic Hilbert space (which we use throughout the remainder of this paper), Schack and 
one of us [Tij introduced a class of quantum baker's maps for dimensions D = 2 N . For these cases, we can model our 
Hilbert space as the the space of N qubits, and the position states can be defined as product states for the qubits in 
the standard basis, i.e., 



kj) = \xi) ® \x 2 ) <8> • • • <8> \x N ) , (7) 



where j has the binary expansion 



N 

j = x l ...x N -G = Y J Xi2 N - 1 (8) 

and qj = (j + 1/2)/ D = 0-x% . . . x^l. 

The connection with the classical baker's map derives from the symbolic dynamics. The bi-infinite strings J2J that 
specify points in the unit square are replaced by sets of orthogonal quantum states created through the use of a partial 
Fourier transform 

G n EE l 2 n ® F 2 N-n n = 0, . . . , N, (9) 

where I2™ is the unit operator on the first n qubits and F 2 N-n is the Fourier transform on the remaining qubits. The 
partial Fourier transform thus transforms the N — n least significant qubits of a position state, 

G n \xi) 8 ■ ■ ■ ® \x n ) ® |m) ® ■ ■ ■ ® \a N - n ) = \xi) ® ■ ■ ■ ® \x n ) ® ^ — V \x n+1 ) <E> ■ ■ ■ <E> \x N )e 27Tlax/2N ~" , (10) 

v2 1 — 

X„ + 1,...,XJV 

where a and x are defined through the binary representations a = 01 . . . q.at_„-1 and x = x n+ \ . . . x^-l. In the limiting 
cases, we have Go = Fjj and Gn = *1- The analogy to the classical case is made clear by introducing the following 
notation for the partially transformed states: 

I aAr_ n . . . a% • xt . . . x n ) ee G n \x x ) ® • • • ® \x n ) ® \ai) ® • • • <g> |aAr_„) . (11) 

For each value of n, these states form an orthonormal basis and are localized in both position and momentum. The 
state \dN-n ■ ■ ■ o>x • x\ . . .x n ) is strictly localized in a position region of width 1/2™ centered at 0-x± . . .x n l and is 
roughly localized in a momentum region of width l/2 N ~ n centered at 0-ai . . . ajv-nl- In the notation of Eq. it 
is localized at the phase-space point la^-n ■ ■ ■ a\ • x\ . . . x n l. Notice that |ojy ■ ■ ■ oi*) = Go |oi) (8 • • • (8 |ojv) is a 
momentum eigenstate and that |» x% . . . xn) = Gn \xx) ® • • • ® \xn) — i\xi) (8> • • • <8> \xn) is a position eigenstate, the 
i being a consequence of the anti-periodic boundary conditions. 

Using this notation, a quantum baker's map on qubits is defined for each value of n = 1, . . . , N by the single- 
iteration unitary operator [T(| 

BN,n = G n _ 1 o S n ° G^ 1 = I Q>N—n ■ ■ ■ &lXx * X 2 . . . X n ) (dN — n ■ • ■ Qi\ * X\X 2 . . . X n \ j (12) 

a:i,...,a:„ ai,...,a^v — ^ 

where the shift operator 5„ acts only on the first n qubits, i.e., i§ n |a;i) ® lajg) <8> • •• <8> |x„) ® |a;„+i) ® • •• <8> |a;jv) = 
| £2) <E> ■ • • ® \x n ) ® |iCi) 8> ) <8 •• • ® |a;jv). Notice that since commutes with G" 1 , we can put Bn, u in the form 

(U^F^^oSn- (13) 

Since Si is the unit operator, it is clear that Bn,i is the Balazs-Voros- Saraceno quantum baker's map ©. 
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(a) « (b) « (c) « 

FIG. 1: Husimi function for each partially Fourier transformed state (II 1 I t when N = 2: (a) n = 2, (b) n = 1, and (c) n = 0. 
The action of the quantum baker's map £2,2 is to map the four states in (a) to the four states in (b), as shown by the numbers 
labeling the states. Similarly, the action of £>2,i is to map the states in (b) to the states in (c). 



We can also write 

B N ,n = hn-i ® fljv-n+1,1 S n , (14) 

which shows that the action of Bn.u is a shift of the n leftmost qubits followed by application of the Balazs-Voros- 
Saraceno baker's map to the N — n + 1 rightmost qubits. At each iteration, the shift map S n does two things: it 
shifts the n-th qubit, the most significant qubit that was subject to the previous application of 5jv- n +i i> out of the 
region subject to the next application of -B/v-n+i,ii and it shifts the most significant (first) qubit into the region of 
subsequent application of -Bjv-n+1,1- 

The quantum baker's map Bm,u takes a state localized at laN-n---ai • xi...x n l to a state localized at 
lajsr—n . . . a\X\ • xi . . . x n l. The decrease in the number of position bits and increase in momentum bits enforces 
a stretching and squeezing of phase space in a manner resembling the classical baker's map. In Fig.^a), (b), and (c), 
we plot the Husimi function (defined as in f° r the partially Fourier transformed states when N = 2, and 
n = 2, 1 and 0, respectively. The quantum baker's map is a one-to-one mapping of one basis to another, as shown in 
the figure. 

One useful representation of our quantum baker's maps, introduced in |10| . starts from using standard techniques 
[2l| to write the partially transformed states <|XU|) as product states: 




_|_ e 27ri(Oa JV _fe + i...a JV -„l)|-|A j ^ 



These input states are mapped by to output states 

I ajv _„ . . . a lXl . x 2 . . . Xn ) = e **(°«i«i-«"-i) 

Y| ) + e 2 " (0aN -' t + 1 - aN -" 1) |l)) ) (g) 4=fl°) +e 2 " (0xiai - ajv -" 1) |l)) 




A 1 f, 



(16) 



These forms show that the quantum baker's map Sjv,„ shifts the states of all the qubits to the left, except the state 
of the leftmost, most significant qubit. The state \x\) of the leftmost qubit can be thought of as being shifted to the 
rightmost qubit, where it suffers a controlled phase change that is determined by the state parameters a\, . . . , apj^ n 
of the original "momentum qubits." The quantum baker's map can thus be written as a shift map on a finite string of 
qubits, followed by a controlled phase change on the least significant qubit. Soklakov and Schack ^lj have developed 
this shift representation into a useful tool. Using an approach based on coarse graining in this representation, they 
investigated the classical limit of the quantum baker's maps. 
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Another useful representation of our quantum baker's maps is the qubit (position) representation 

Bn,ti = 2 N-n+l E E \ X ^){ X A <8> ••• ® \x n ){x n -l\ ® ]^l><a;„l ® \z 2 ){yi\ ® ••• ® |«JV-n+l)(yjV-n| 



x exp 



)JV- 



-((j + 1/2) (Z + 1/2) + 2 JV -"x 1 (Z + 1/2) - 2(j + 1/2) (fc + 1/2)) 



(17) 



where 



7V-n JV-n N—n+1 

fe=i fc=i fe=i 

Soklakov and Schack |Tl) have found simplified forms of this qubit representation, suitable for asymptotic analysis of 
the classical limit. 

The classical limit for the above quantum baker's maps has also been investigated in using an analysis based 
on the limiting behavior of the coherent-state propagator of Bn.u- When D = 2 N — > oo, the total number of qubits 
N necessarily becomes infinite, but one has a considerable choice in how to take this limit. For example, we could 
use only one position bit, thus fixing n = 1, and let the number of momentum bits N — 1 become large. This is 
the limiting case of the Balazs-Voros-Saraceno quantization. There is, however, a wide variety of other scenarios to 
consider, e.g., n = N/2 or n — 2N/3 — 1 as N — ► oo. In it was shown that provided the number of momentum 
bits N — n approaches infinity, the correct classical behavior is recovered in the limit. If the number of momentum 
bits remains constant, e.g., n = N as N — > oo, a stochastic variant of the classical baker's map is found. These results 
are consistent with those obtained previously by Soklakov and Schack . 

The special case n = N, which does not limit to the classical baker's map, has other interesting properties. Although 
all finite-dimensional unitary operators are quasi-periodic, the quantum baker's map -Bjv.jv is strictly periodic, 

B NJV = 1 j ( 18 ) 

as we show below. All its eigenvalues, therefore, are 4iV-th roots of unity i.e., of the form e 7r4fc / 2Ar ; and, hence, 
there are degeneracies when N > 4. This represents a strong deviation from the predictions of random matrix theory 
[2^,|23,|2^|. The property (|T5|l can be easily shown by noting that Bjv.iv = — iGn-i SW = — i(i 2 «-i <X> F%) o §n; i.e., 
Bn.n is a shift followed by application of the unitary 

U = -iF 2 = -L( e -W4(|o)( | + |l)(l|) + e"/ 4 (\0)(l\ + |1)(0|)) = e -WVM-/4) ( 19) 

to the least significant qubit. On product states, the action of Bn,n can be written explicitly as 

B n ,n\$i) ® IV>2> ® ■ ■ ■ <8> \ip N ) = ^2) <E> ... |Vat> <E> U\i/ji) . (20) 

Since = 1, we get the property JTSJ. One can now also see that Bn.n cannot entangle initial product states. 

The eigenstates of U are (|0) + |l))/\/2~, with eigenvalue +1, and (|0) — |l))/\/2, with eigenvalue —i. For the 
discussion in this paragraph, label these eigenstates by their eigenvalues a = We can construct eigenstates of 

Bn,n from tensor products of these eigenstates. Let P denote the period of a string a± . . . aw under cycling; for the 
the corresponding product state \ip) — \a±) (8 • • • ® |ajv), P is the smallest positive integer such that S^\tjj) — \ip). It 
is now easy to show that the eigenstates of Bn.n have the form 



1 P ~ 1 

*) = -jp t" 1 " ' a p\p ka ^ ■ ■ ■ a kS k N \ai) ® ■ ■ • ® |ajv> , (21) 



k=0 



where [a)p denotes a P-th root of a. The eigenvalue of the state l|21|) is [a\ ■ ■ -ctp]p. Notice that the product state 
|a) (8Ar (a — 1, —i) is always an eigenstate, with eigenvalue a. As an example, the eigenstates of i?2,2 are 

|1>®|1>, |-i)®|-<), -^=(|l)®|-i)+ e "/ 4 |-i)®|l)), -^(ll^l-O + e- 3 ™/ 4 !-,)®!!)) , (22) 
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with eigenvalues 1, —i, e -7 ™/ 4 , and e 37 ™/ 4 , respectively. 

When n < N, the action of the quantum baker's map is similar to l|2(J|l . but with a crucial difference. After the qubit 
string is cycled, instead of applying a unitary to the rightmost, least significant qubit, a joint unitary is applied to all 
of the the N — n + 1 rightmost qubits. As discussed above, this joint unitary can be realized as controlled phase change 
of the least significant qubit, where the control is by the state parameters eti,. . . ,aN- n of the original momentum 
qubits. This controlled phase change means that initial product states become entangled. The resulting entanglement 
production is the subject of this paper and is investigated in Sec. IIVI Since the entangling controlled-phase change 
involves an increasing number of qubits as n decreases from n = N to n = 1 (the Balazs-Voros-Saraceno map), we 
might expect the entanglement to increase as n ranges from TV to 1. What we find, however, is that all the maps for 
n not too close to N are efficient entanglement generators, but that the greatest entanglement is produced when n is 
roughly midway between TV and 1. 

To calibrate our entanglement production results, in the next section we establish how much entanglement to expect 
for pure states chosen randomly from the Hilbert space. With the entanglement of typical states quantified, we have 
a standard against which to compare the entanglement produced by the quantum baker's map. We might expect 
the quantum baker's maps to be good at creating typical states in Hilbert space, and our work on entanglement 
production can also be regarded as a way of investigating this expectation. 

III. ENTANGLEMENT OF TYPICAL PURE STATES 

Quantifying the amount of entanglement between quantum systems is a recent pursuit that has attracted a diverse 
range of researchers |2jJ, 1211 liil Li2| ■ The best understood case, not surprisingly, is the simplest. It is generally accepted 
that when a bipartite quantum system is in an overall pure state, there is an essentially unique resource-based measure 
of entanglement between the two subsystems. This measure is given by the von Neumann entropy of the marginal 
density operators |33l l34j . Thus an investigation into typical values expected for the entropy of entanglement seems 
a worthwhile endeavor in its own right 0, |3l|, which we undertake in Sec. IIII Al Notice that given an A-qubit 
quantum-baked state, there are 2 — 1 different possible partitions into the two subsystems, and hence, 2 — 1 
different entropies of entanglement to consider. 

Another well understood case is the pairwise entanglement of two qubits in an overall mixed state. When a bipartite 
quantum system is in a mixed sta te, prop osals for measuring the ent ang lement include the entanglement of formation 
pfl l33| , distillable entanglement j.'i'ii |30[ | , and relative entropy [3?l l38j . For pure states all these reduce to the von 
Neumann entropy, but the first is the best understood for mixed states. In the case of two qubits in a mixed state, 
an exact expression for the entanglement of formation exists in terms of another measure called the concurrence 
[30l 13^. |40| . Thus another entanglement measure to consider could be the concurrence that results after all but two 
qubits are traced out of an A-qubit quantum-baked state. Section [ill Bl is concerned with this pairwise entanglement. 

Unlike the above special cases, quantifying the amount of multipartite entanglement in a general multipartite state 
remains far from being completely understood. There have, nevertheless, been a number of proposals for such a 
measure. In Section IIII CI we investigate two of these. We must stress, however, that no single measure alone is 
enough to quantify the entanglement in a multipartite system. As the number of subsystems increases, so too does 
the number of independent entanglement measures. 

A. Bipartite pure-state entanglement 

Consider a bipartite quantum system with Hilbert space Ha <8> T~Cb of dimension pv, where \x and v are the 
dimensions of subsystems A and B, with fi < v. A joint pure state p = |V')( 7 /'I nas a Schmidt decomposition 
\ip) = X)f=i \/Pi\ a i) ® |ki)> where |a.;) and \bi) are orthonormal bases spanning Ha and Hb- If we sample random 
pure states according to the unitarily invariant Haar measure, then the Schmidt coefficients < pi < 1 obey the 
distribution |4l| 

II (Pi-ptfliPr^Pk, (23) 

1=1 / l<i<j<n k=l 

where A is a normalization constant. 

Considered as eigenvalues of the marginal density matrices pa = tr^p and ps — tr^p, the Schmidt coefficients 
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give the von Neumann entropy S of each subsystem 

M 

S = Sa = S B = - tr p B lnps = - J^P; lnpj. (24) 

i=i 

As remarked above, the von Neumann entropy is generally considered to be the unique resource-based measure of 
entanglement for a bipartite quantum system in an overall pure state. Given the distribution (|23[) . an expression for 
the average entropy can be calculated: 

S,^(S)= £ l-^. (25) 

k=v+l 



This succinct formula was conjectured by Page |42| and later proved by others |43| |44j, |45j (see also |19ll35l l46l |47 | ). 

An expression for the average purity, 

R = R A = Rb = trp| = Y^Pf . ( 26 ) 
had been calculated much earlier by Lubkin ^| : 

= (R) = . (27) 

liv + 1 

The purity provides the first nontrivial term in a Taylor series expansion of the von Neumann entropy about its 
maximum and because of its simplicity, is much easier to investigate analytically. For these reasons we restrict our 
attention to it. 

One can also define a linear entropy in terms of the purity: 

S L = (3(1 - R) . (28) 
We choose the normalization factor so that < Sl < 1, i.e., (3 = /i/(/i — 1). The average linear entropy is 

which shows that for any division into two subsystems, when the overall dimension [iv is large, a typical state has 
nearly maximal entanglement. 

Ideally, we would like an expression for the complete probability distribution P(R) dR of the purity for random pure 
states. This function cannot, in general, be calculated analytically, so we are forced to settle for formulae describing 
a few of the cumulants. For subsystems of even moderately high dimension, the cumulants we calculate are sufficient 
to describe accurately the entire distribution P(R). In deriving these cumulants, we follow the work of Sen |4q . 

Consider the second moment 

(r 2 ) = [ ir P ? P fp(p)dp, (30) 

where p = (pi, . . . ,p M ) and dp = dpi . . . dp^. We first remove the obstacle of integrating over the probability simplex 
by noting that 

t 1 

Q(q)dq= [] (q i -q j ) 2 l[e-'"'q v k - fi dq k =Ne- r r^- 1 P(p)dpdr, (31) 

l<i<j<fj, k=l 

where the new variables qi = rpi take on values independently in the range [0,oo). Integrating over all the values of 
the new variables, we find that the normalization constant is given by N = Q /T(fiis), where Q = J Q(q)<iq. Similarly, 
we find that 



(32) 
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thus determining the desired second moment in terms of integrals over Q(q). 

Now notice that the first product in Eq. (|31|) is the square of the Van der Monde determinant 



A( q ) = J] = 



1 . 


. 1 








<7i • 


■ 9m 




rf(ffi) ■ 




u-l 

ll ■ 


. q^ 1 




r?_i(9i) ■ 





(33) 



The second determinant in Eq. (|33[1 follows from the basic property of invariance after adding a multiple of one row to 
another, with a = v — fx and the polynomials r?(q) = k\L^(q) judiciously chosen to be rescaled Laguerre polynomials 
|49| . satisfying the recursion relation 



rUl) = rt + \q) - kr^l(q) = £(-l)< J )k(k - 1) . . . (k - i + l)r^(q) 

and having the orthogonality property 

dq e~iq a rl (q)rf (q) = T(k + l)T(a + k + l)6 kl . 
These facts in hand, we can evaluate 

Q = / A(q) 2 Y[e~^q%dq k 
J k=i 

E <<••••< <,.•••., II / d?*e- ft gg'rf fc _ 1 (g fc )rSL_i(«*) 



fc = l ' 



= e 4...vii r fe) r ("+u) 
ti,...,^ fc=i 

= /*! ]J r(fc)r(a + fc) . 



fc=i 

Elaborations of this calculation lead to the following formulae: 

n-i 



E/^Q(q)rfq = ^Ef(fcTT 



4 4 » 



E / = Q E 



r(fc + i)r(a + fc + 1) ' 

fcl^o r ( fc + ^ r ( a + k + + 1 ) r ^ + 1 + ^ 



fe=0 



Here 



r(fc + i) E (~i) 4+r 

i,r— 



3 \ 3 



T(a + j + l-r+ 1)1(1 -l)...(l-r+ !)*,_,.,*_< 



(34) 



(35) 



(36) 



(37) 
(38) 



(39) 



where the final form follows from Eqs. I|34|l and (|35|l . Evaluating the sums in Eqs. (|37|) and l)38|) leads to the 
simplification 

E / 9?9jQ(q) ^ = + ^)(^ 2 + + 5^*1/ + 5) + O - l)(u - 1)0 + v - + v - 2)] . (40) 



9 




0.94 0.95 0.96 0.97 0.98 0.91 0.92 0.93 0.94 0.95 

S r ( 5 = 0.9650 ± 0.0062 ) S. ( 5= 0.9339 ± 0.0058 ) 

(c) (d) 



FIG. 2: Histograms for the probability distribution P(Sl), numerically calculated from 1 million random states for the several 
subsystem dimensions n and v possible in a 256-dimensional Hilbert space. The exact means (vertical solid line) and standard 
deviations (vertical dashed line), along with Gaussian (dotted curve) and Airy function (dash-dotted curve) approximations to 
the distribution are also plotted. For /i = 2, the exact distribution 1491 is plotted as the solid curve. Notice that for /i not too 
close to 2, the Gaussian approximation is sufficient and, as noted in the text, the linear entropy of a typical state, though close 
to maximal, is localized away from maximal. 



Using Eqs. (|30|) and l|32l) . we now obtain 



2 = Qi + v)i^ + v 2 + hjiv + 5) + {ji - \){y - l)(n + v- l)(/i + v-2) 
[ ' (fu/ + 3)0*1/ + 2)(jh/ + 1) ' 1 ' 



The variance is then given by 



ms{ ^- {Rf = T -^)i^- T? . ,42, 



Using the same methods, one can also derive an expression for the third cumulant. Due to the complexity of these 
calculations, we only state the final result: 

«i? 3 )) (R*) - 3(i?)(i? 2 ) + 2(Rf = - 8{ £' 1K ^ 1 ^ 1/)0 Z7 5) rTs ■ ( 4 3) 
vv " v ' v /v ' w (/i^ + 5)(/ii/ + 4)(/ii/ + 3)(/ii/ + 2)(/ii/+ l) 3 v ' 

Translating our results to the linear entropy, we have 

iS L )) = (5 £ )^/3 (/i ~ 1) ^~ 1) (=o), (44) 
fiv + 1 
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N=3 c = 0.369 ±0.176 



FIG. 3: Pairwise mixed-state entanglement in random multi-qubit pure states with a total of 3, 4, 6, and 8 qubits. The plots 
show the probability distribution of the quantity c(p) = Ai — A2 — A3 — A4. The concurrence is C(p) = max{0, c(p)}. The means 
and standard deviations of the distributions are also given. The approximate probability to encounter pairwise entanglement 
between a particular pair of qubits in states with 3, 4, 6, and 8 total qubits is 1, 0.76, 0.006 and 0, respectively. The distributions 
were numerically calculated using 1 million random states. 



((St)) 



P 2 ((R 2 )) 
-P 3 ((R 3 )) 



(= b) 

(= c) 



(45) 
(46) 



These can be used in an approximation to the cumulant generating function and, hence, to the probability distribution 
itself: 

1 f°° 

P(S L )dS L ~ 7- \ duj exp [-iS L w + aiu + b{iu) 2 /2\ + c{iuf /3!] dS L 

z ™ J — oo 

, b(S L -a) 



= |2/c| 1/3 exp 



Ai 



(2/c 



,1/3 



Sl — a 



dSh ■ 



3c 2 c 

When the overall dimension \xv is large, the standard deviation of the linear entropy is given approximately by 

5 



((si)) 



2\\l/2 



[iv + 1 V M 



1 



1 



1 



(47) 



(48) 



Comparing this with the average linear entropy shows that the bipartite entanglement of a typical pure state, though 
close to maximal, is nonetheless localized away from maximal as long as ((S^)) 1 ^ 2 is somewhat smaller than (// + 
l)/(fj,v + 1), i.e., as long as fj, is somewhat larger than 2. 

In Fig. |3we display numerical calculations of P(Sl) for the several ways of dividing a 256-dimensional Hilbert space 
into two subsystems. These numerical calculations used 1 million random states. We also plot the means (vertical 
solid line), standard deviations (vertical dashed line), the Airy function approximation (|47|l (dash-dotted curve), and 
the Gaussian approximation (dotted curve). For the special case \x = 2, the exact probability distribution is drawn 
(solid curve): 



P(S L ) dS L = 2T Xt 1/2 1 V^SlS L "- 2 dS L (ji = 2) 
^ttT{v - 1) 



(49) 



Note that the distributions are highly localized and that for fi somewhat larger than 2 in a high-dimensional overall 
space, the Gaussian approximation is sufficient. 



B. Pairwise mixed-state entanglement 

A numerical study of pairwise (mixed-state) entanglement in multi-qubit systems has already been published |50j | , 
making the results in this section somewhat redundant. Our choice for an entanglement measure is the concurrence 
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[30| of a two-qubit density operator p, given by 

C(p) = max{0, A x - A 2 - A 3 - A 4 } (50) 

where Ai > A2 > A3 > A4 are the square roots of the eigenvalues of p(a y a y )p*{a y £g> <j y ). The complex conjugation 
is taken in the standard qubit basis. The concurrence takes values in the range [0, 1], with a pair of qubits being 
entangled if and only if C(p) > 0. 

The concurrence provides an explicit formula for the entanglement of formation 

E f (p) =in{ J^PjS^j), (51) 
j 

where the infimum is taken over all pure-state decompositions p — J2j Pj\^ J j}(' l l } j\> an d S(tpj) is the subsystem von 
Neumann entropy of the bipartite pure state ipj . In the case of two qubits, 

E f (p)=£(C(p)) , (52) 

where £ is defined in terms of the binary entropy function h(x) = — xlogx — (1 — x) log(l — x) by 

e(C ,_*(i±4E£). ,53) 

The function 8(C) is monotonically increasing for < C < 1, and hence the concurrence is a good measure of 
entanglement in its own right. 

To apply the concurrence as a measure of pairwise entanglement in iV-qubit pure states, we first trace out N — 2 
of the qubits and then use the above formulae on the resulting two-qubit density operator. In Fig. [2 we have plotted 
the probability distribution for the quantity 

c(p) = Ai - A 2 - A 3 - A 4 (-1/2 < c < 1) , (54) 

when the iV-qubit pure states are sampled from the Haar distribution for N — 3, 4, 6, and 8. Notice that the 
probability of finding pairwise entanglement between any pair of qubits decreases rapidly as N increases. In contrast, 
the preceding subsection shows that as N increases, the bipartite entanglement of a typical state is close to maximal 
no matter how the overall system is sliced into two parts. Taken together, these results mean that the entanglement 
in a typical state of many qubits is mainly multipartite entanglement shared among many of the qubits, rather than 
pairwise entanglement. 



C. Multipartite entanglement 

We now investigate two proposals for a measure of multipartite entanglement, the measure Q of Meyer and Wallach 
|5l| and the n-tangle of Wong and Christensen |52l |. In general, as the number of subsystems increases, an exponential 
number of independent measures is needed to quantify fully the amount entanglement in the multipartite system. 
Consequently, neither of the following entanglement measures can be thought of as unique. Different measures capture 
different aspects of multipartite entanglement. 

The Meyer- Wallach measure, Q(ip), which can only be applied to multi-qubit pure states, is defined as follows. For 
each a = 1,...,N and 6 £ {0, 1}, we define the linear map i a (b) : (C 2 )® N —> (C 2 )®^" 1 through its action on the 
product basis, 

i a (b)\xx) (g) . . . (g) \x N ) = 5 bXa \xi) ® . . . (g> \x a -i) <g> \x a+ i) (g) . . . <g) |a;jv) • (55) 

Next let 

D(il>,4>) =X)l^-^ii 2 . (56) 

be the square of the wedge product of two vectors and |</>), where the ipj are the coefficients of the state \tp) in 
the product basis, 

2 m rn 

\tp) =^2i>j\xi>® \x m ) , j = x 1 ...x m -0 = J2 x i 2m ~ l ■ ( 57 ) 
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FIG. 4: Multipartite entanglement in random multi-qubit pure states with a total of 2, 4, 6, and 8 qubits. The probability 
distributions are for the multipartite entanglement measures Q and tjv. For comparison, the distributions for the subsystem 
linear entropy Sl are also plotted in the case of equal subsystem dimensions (ft = v = yD). The analytic values of the means 
and standard deviations are also shown. The distributions are based on numerical calculations using 1 million random states. 



The Meyer- Wallach entanglement measure is then 



iV 



(58) 



a=l 



Meyer and Wallach have shown that Q is invariant under local unitary transformations and that < Q < 1, with 
Q(ip) = if and only if \ijj) is a product state. It was recently shown by Brennen |5^| that Q is simply the average 
subsystem linear entropy of the constituent qubits: 



(59) 



Here pk is the density operator for the k-ih qubit after tracing out the rest. Hence we should expect this measure to 
agree qualitatively with the subsystem entropies. 



13 




(a) 






50 


100 


150 


200 








iterations 








(b) 



(c) 




200 




50 100 150 200 

iterations 




(d) 



100 150 

iterations 



200 



0.95 




0.85 



(e) 



100 

iterations 



200 




(0 



100 

iterations 



200 



FIG. 5: Dynamical behavior of the subsystem entropy, (a) and (b) plot the von Neumann entropy and the linear entropy, 
respectively, when N = 8, the partition divides the four rightmost qubits from the four leftmost, and the initial state is 
j 00000000). (c) is the same as (b), except that the partition divides the single rightmost qubit from the others, (d) and (e) are 
also the same as (b), but with initial state |001001 11) in (d) and maximally entangled initial state (1651 in (e). (f) shows the 
case of TV = 4 qubits, a partition that separates the two leftmost qubits from the two rightmost, and an initial state jOOOO). 
The different quantum baker's maps are colored according to the key in (d). The dashed lines show the average entanglement 
predicted for random states. 
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Using, for example, the Hurwitz parametrization |54l l55j| of a random unit vector in C D , where the dimension 
D = 2 N , one can calculate the mean and variance of Q for random states: 

(Q) = §^f , (60) 
((Q 2 )) = (Q 2 ) (Q) 2 = 6( ° - 4) + — (61) 

The mean was also calculated independently in |56j |. and given the relationship (|59|) . it can be easily checked using 
Lubkin's formula for the average purity (|27p. For large N, the mean and standard deviation are given approximately 
by (Q) w 1 — 3/1? and ((Q 2 )) 1 ^ 2 ~ (1/ 'D)y/6/N , indicating that the Meyer- Wallach entanglement of a typical state is 
very nearly maximal, but suggesting that this entanglement is weakly localized just below maximal. 

In the case of two qubits, the square of the concurrence is often referred to as the tangle. A generalization of the 
tangle to three qubits was defined by Coffman et al. 57]. Wong and Christensen proposed another generalization 
valid for an arbitrary even number of qubits. For pure states of N (even) qubits, it is defined as 

T N ^) = Maf N \r)? (iVeven). (62) 

Wong and Christensen were able to show that tjv is an entanglement monotone and also to generalize its definition 
to mixed states. From the above definition, it is easy to see that < t n < 1. A cat state, (jO)®^ + \1)® N ) / y/2, has 
maximal entanglement by this measure. 

One can calculate the mean and variance of tn for random states: 

(63) 



= (/; + 3) (/;'■■■ ■ m > 



For large N, both the mean and the standard deviation of tn are given approximately by 2/D, indicating that a 
typical state does not have whatever sort of multipartite entanglement is characterized by t/v ■ 

In Fig. 0] we plot the probability distributions of Q and tjv after sampling 1 million random multi-qubit pure 
states with a total of 2, 4, 6, and 8 qubits. Also plotted, for comparison, are the corresponding distributions for 
the subsystem linear entropy Sl, with the two subsystem dimensions chosen to be equal. When N = 2, the three 
different measures are equivalent, each having the distribution P(E)dE = (3/2)vl — E dE (E = Sl,Q,tn)- Notice, 
however, that the behaviors of the two multipartite measures diverge dramatically as we increase the total number of 
qubits. According to the Meyer- Wallach measure, multipartite entanglement increases as we increase the number of 
qubits. This agrees with the bipartite measure Sl, the only difference being that the Meyer- Wallach entanglement 
of a typical state is closer to maximal than the linear entropy as N gets large. In contrast, the n-tangle of a typical 
state decreases as N increases. As noted above, tjv seems to describe some sort of multipartite entanglement that 
becomes rare as the number of qubits increases. 



IV. ENTANGLEMENT PRODUCTION IN THE QUANTUM BAKER'S MAPS 

We now report the results of a detailed numerical study of entanglement production for the quantum baker's maps. 
All of the results reported in this section are for the 8-qubit baker's maps (except for Fig.|SJf), which displays results 
for N = 4 qubits). In the figures, the maps for different values of n are distinguished by a color coding given in 
Fig- Eld). The extreme map -B/v.at is not included in the plots because it does not produce any entanglement. 

First consider the bipartite measures of entanglement. In Fig. we plot the dynamical behavior of these measures 
for several interesting initial states. Figure [5fa) shows the behavior of the subsystem von Neumann entropy S when 
N = 8 and the bipartite division is between the four least significant qubits (rightmost) and the four most significant 
(leftmost); in terms of the phase-space description of the baker's maps, we are considering entanglement between the 
fine and coarse position scales on phase space. The initial state in Fig. [Sfa) is the product state 1 00000000). Note 
the ranking among the different quantum baker's maps. The maps -Bs,4 (red) and Bs,3 (dark blue) achieve saturation 
values closest to that for typical states (dashed line). The original Balazs-Voros-Saraceno quantization has n = 1 
(black). Figure |5jb) is the same as (a), except now the linear entropy Sl is our entanglement measure (as it remains 
in the remaining parts of Fig. [SJ. The two different definitions of entropy show the same qualitative behavior. In 
Fig. [5fc), we again use linear entropy, but switch to a bipartite division that separates the single rightmost qubit 
from the rest. The different saturation levels attained by the different quantum baker's maps, though less intelligible, 
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FIG. 6: Distributions for the subsystem linear entropy, relative to a partition that divides the four least significant qubits from 
the four most significant, after 5, 20, and 500 iterations of the quantum baker's maps starting with a uniform distribution of 
16 000 product states (delta function centered at zero). A total of N = 8 qubits is used. The dashed-dotted lines are the means 
of the distributions, and the mean and standard deviation for random states are the solid and dashed lines, respectively. 
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FIG. 7: Evolution of the means and standard deviations for the subsystem linear entropy under the same conditions as in 
Fig. HJ The dashed lines are the mean and standard deviation for random states. The middle, magnified plot displays the 
ranking of the baker's maps in terms of entangling power, i.e., 4,5,3,2,1,6,7. Error bars (2A5l/V16 000) are included in the 
middle plot to indicate the expected sampling error for the means at saturation; these error bars show that the rankings in 
terms of entangling power are real, not statistical artifacts. 



are still discernible. If we return to our original partition and change the initial state to |00100111), however, the 
differences in the saturation levels disappear almost altogether, as is shown in Fig. Efd). The quantum baker's map 
with n = 7 stands out in all cases, saturating at a value well below the other maps; we should recall that the map 
£?8, 7 is closest to the map Bs,8, which does not entangle at all. In Fig. GJe), again using the original partition, we 
start in the entangled state 

j ^2 \x\X2XzXiXxX2XzX4) , (65) 

Xi,X 2 ,X3,X4 

which is maximally entangled with respect to the original partition. In this case the maximal initial level of entangle- 
ment is destroyed by the quantum baker's maps. Figure^Jf) shows the case with N = 4 qubits, a partition separating 
the two leftmost qubits from the two rightmost, and an initial state of 1 0000) . 

In view of the variety of behaviors exhibited by subtle differences in the initial states [compare Fig.|5fb) toOJd)], if 
we are to study the intrinsic entangling properties of the quantum baker's maps and not properties conditioned on a 
particular initial state, we need an approach that treats all initial states of a certain type on the same footing. Such 
neutrality can be achieved by takin g an average over, for example, the set of all product states. This approach was 
used to define the entangling power [58( of a unitary operator, and we adhere to it for the remainder of this section. 
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FIG. 8: Evolution of the means and standard deviations for the quantity c of Eq. 1541 . where all but the single leftmost and 
rightmost qubits are traced out. The concurrence for the remaining two qubits is C = max{0, c}. The averages are taken with 
respect to 16 000 uniformly distributed initial product states for N — 8 qubits. The dashed lines are the numerically calculated 
mean and standard deviation for random states. Error bars (2Ac/vl6 000) are included in the middle plot to indicate the 
expected sampling error for the means at saturation. 



Consider starting with a uniform distribution of product pure states, each member being a tensor product of N = 8 
random single-qubit states. We now "bake" entanglement with the seven quantum baker's maps. The remaining 
figures plot the amount of baked entanglement as quantified by the several entanglement measures discussed in 

sec, urn 

Figure El plots the distributions of the subsystem linear entropy, relative to a partition that divides the four least 
significant qubits from the four most significant, after 5, 20, and 500 iterations of the quantum baker's maps. All the 
distributions start in a delta function centered at zero entanglement, but quickly spread while advancing toward the 
value predicted for typical states. The mean and standard deviation for random states are the solid and dashed lines, 
respectively. At 500 iterations, well after saturation (~ 100), the mean entanglement for our quantum-baked states 
(dash-dotted lines) fall short of the value predicted for random states. Notice that the variances are also still quite 
large. For clarity, the means and standard deviations are plotted separately in Fig. [7| The ranking of the different 
quantum baker's maps in terms of entangling power — 4,5,3,2,1,6,7 — is now only apparent after magnification at the 
saturation level (middle plot). It was found that this ranking is preserved when the partition is changed. If the 
total number of qubits is varied, similar behavior also results, and we conclude that the quantum baker's maps are, 
in general, good generators of bipartite entanglement, with the highest levels produced when n is roughly midway 
between N and 1. A total of 16 000 initial states were used in this simulation. The same number is used for those 
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FIG. 9: Evolution of the means and standard deviations for the Meyer- Wallach multipartite entanglement measure Q. The 
averages are taken with respect to 16 000 uniformly distributed initial product states for N = 8 qubits. The dashed lines are 
the numerically calculated mean and standard deviation for random states. Error bars (2AQ/V16 000) are included in the 
middle plot to indicate the expected sampling error for the means at saturation. 

that follow. 

We next consider the pairwise entanglement as given by the concurrence, or more specifically, by the quantity c(/3) 
of Eq. H54JI. We only display results for the case where all but the single leftmost and rightmost qubits are traced out. 
The other cases are similar. Again, starting in a uniform distribution of 16 000 product states (delta function centered 
at c = 0), we bake entanglement into our states. In this case, however, the pairwise entanglement is not maintained. 
After a short rise, c(/5) quickly falls to negative values, retreating toward the mean numerically calculated for random 
states, as shown in Fig. EI Not surprisingly, the quantum baker's maps respect the rarity of pairwise entanglement in 
multi-qubit states. The ranking among the different quantum baker's maps is also preserved. 

In Figs. 151 and llOl we plot the corresponding evolutions for multipartite entanglement measures, the Meyer- Wallach 
measure Q and Wong-Christensen tjv- The means and standard deviations of Q behave similarly to the bipartite 
measures of entanglement. The ranking among the different quantum baker's maps is again preserved. This is not 
surprising given the relationship l|59|l . In the case of tjv, however, it is difficult to discern any useful information, 
presumably because this measure only describes entanglement of a very special type, e.g., that in ./V-qubit cat states. 
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FIG. 10: Evolution of the means and standard deviations for the Wong-Christensen tangle rjv. The averages are taken with 
respect to 16 000 uniformly distributed initial product states for N — 8 qubits. The dashed lines are the numerically calculated 
mean and standard deviation for random states. Error bars (2Atjv /V 16 000) are included in the middle plot to indicate the 
expected sampling error for the means at saturation. 



V. DISCUSSION AND CONCLUSION 



The numerical calculations of the previous section show that the quantum baker's maps are, in general, good at 
creating multipartite entanglement amongst the qubits. It was found, however, that some quantum baker's maps 
can, on average, entangle better than others, and that all quantum baker's maps fall somewhat short of generating 
the levels of entanglement expected in random states. This might be related to the fact that spatial symmetries in 
the baker's map allow deviations from the predictions of random matrix theory Q. Such deviations are apparent in 
the statistics of the eigenvectors of Bn,u and might also taint the randomness of our quantum-baked states. In this 
light, an entanglement measure might, in fact, provide a reasonable test for the randomness of states produced by a 
quantum map. 

In our case, entanglement amongst the qubits relates directly to correlations between the coarse and fine scales of 
classical phase space. Although we have only considered entanglement in position, the Fourier transform provides 
a means to investigate entanglement in momentum, and the partial Fourier transform H1(J|) might be used for all 
intermediate possibilities. These qubit bases, while naturally embedded in the construction of the quantum baker's 
maps, might also be applied to other maps of the unit square and, hence, to entanglement production in general for 
quantized mappings of the torus. 

We should expect high levels of entanglement creation in quantum maps that are chaotic in their classical limit. 
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FIG. 11: Saturation value (Q) of the Meyer- Wallach entanglement measure for all quantum baker's maps up to a total of 
N = 10 qubits. This quantity is approximated by taking a time average over the (512fc)-th iterates (k — 1, . . . , 100) and using 
4 000 uniformly distributed initial product states. 

Such maps have a dynamical behavior that produces correlations between the coarse and fine scales of phase space. 
This behavior is described classically in the form of symbolic dynamics. Investigations into entanglement production, 
using the above product bases, allow us to characterize the quantum version of such correlations. To develop a 
complete picture, however, the entangling properties of regular systems must first be addressed. The possibility 
of noncntangling quantum maps, such as Bn,n, which produce stochasticity in their classical limit, must also be 
addressed. 

We can apply our results to a preliminary investigation of the relation between entanglement production and the 
classical limit. As remarked previously, sequences of quantum baker's maps for which the number of momentum 
bits, N — n, is held constant do not approach the classical baker's map in the limit N — > oo, but instead give 
rise to stochastic variants. To relate this behavior to entanglement production, consider the time average of the 
Meyer- Wallach entanglement measure Q, say, 



where, as before, the average ( • ) is taken over a uniform distribution of initial product states. This quantity provides 
the long time saturation value of Q and is plotted for all quantum baker's maps up to total of N = 10 qubits in 
Fig. Note that although one observes a drop in the levels of entanglement for a small number of position bits, n, 
there is no apparent connection between the level of (Q) and the advent of spurious stochastic limits when the number 
of momentum bits, N — n, remains constant as N increases. One might have expected the saturation value to be 
suppressed in such limits, but this appears not to be the case. We tentatively conclude that entanglement production 
in the qubit bases is unrelated to the creation of a stochasticity in the classical limit. A similar picture emerges when 
the subsystem entropy is used as the entanglement measure. 

In conclusion, we have found that the quantum baker's maps are, in general, good at generating multipartite 
entanglement amongst qubits. Given the relation between our qubit bases and classical phase space, this behavior 
should be expected to arise whenever such a quantum map is chaotic in its classical limit. 




(66) 



k=l 
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